Time transformation for random walks in the quenched trap model 
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Diffusion in the quenched trap model is investigated with an approach we call weak subordination 
breaking. We map the problem onto Brownian motion and show that the operational time is 
Sa ~ X]^-cx)(^^)™ where Ux is the visitation number at site x . In the limit of zero temperature we 
recover the renormalization group (RG) solution found by Monthus. Our approach is an alternative 
to RG capable of dealing with any disorder strength. 
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Random walks in disordered systems with diverging 
expected waiting times have attracted vast interest over 
many decades [l[ Q . Two approaches in this field are the 
annealed continuous time random walk (CTRW) model 
and the by far more challenging quenched trap model 
(QTM). Starting in the seventies, the Scher-MontroU 
CTRW approach was used to model sub-diffusive photo- 
currents in amorphous materials Q. Bouchaud showed 
that the trap model is a useful tool for the description of 
aging phenomena in glasses 0, 0] ■ More recently these 
models were used to describe non self averaging Q and 
weak ergodicity breaking 0, 01 which are important for 
the statistical description of blinking quantum dots [1] 
and diffusion of single molecules in living cells 0. 

This manuscript presents a new approach for random 
walks in a fixed random environment. With physical ar- 
guments [H, [13, U I and rigorous mathematics 12 , we 
know that the QTM in dimensions d > 2 is expected 
to qualitatively behave like its corresponding mean field 
CTRW, the latter being exact when d oo. For a ran- 
dom walk in a quenched disordered system intricate cor- 
relations induced by multiple visits to the same site make 
the problem non-trivial and interesting. For that reason 
RG methods Il4l| were used to tackle this problem. 
With RG Machta [10||_found the scaling exponents of the 
QTM and Monthus 14] investigated the diffusion front in 
the limit of zero temperature (see details below). While 
the RG is powerful it has its limitations: a simple ap- 
proach which predicts the diffusion front is still missing. 

We provide the long sought after breakthrough in the 
statistical analysis of sub-diffusion in the QTM. Our ap- 
proach is based on a novel time transformation. It is well 
known that one may decompose the CTRW process into 
ordinary Brownian motion and a Levy process, an ap- 
proach called subordination [15]. In this scheme normal 
Brownian motion takes place in operational time s. The 
disorder is effectively described by a Levy time transfor- 
mation from operational time s to laboratory time t (see 
some details below) . This method is not generally suited 
for random walks in quenched environments since it uses 
the renewal assumption. In a quenched environment this 
very strong assumption implies that a given lattice site 
is visited only once along the path of the random walker. 



i.e. it neglects correlations. So a new approach capable 
of dealing with quenched disorder is now investigated. 
We focus on the one dimensional case since then the de- 
parture from mean field is the strongest. At the end of 
the Letter we explain how to extend our results to other 
interesting cases. 

Quenched trap model [l|, [3, llBl • We consider a ran- 
dom walk on a one dimensional lattice with lattice spac- 
ing equal one. For each lattice site x there is a quenched 
random variable Tx which is the waiting time between 
jump events for a particle situated on x. After waiting 
for a period Tx the particle jumps to one of its two near- 
est neighbors with equal probability. The particle starts 
on the origin x — ai time t — 0, waits for time tq, then 
with probability 1/2 jumps to x = 1 (or x — —1), waits 
there for ti (or r_i) then if the particle returns to x = 
it waits for a time interval tq etc. The {rj;}s are posi- 
tive independent identically distributed random variables 
with a common probability density function (PDF) 
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for Tx oo and < a < 1 . Hence the Laplace transform 
of the waiting time PDF is tp{u) ^ 1 — Au°' + ■ ■ ■ when 
w — >■ 0. As well known [l[ the QTM describes a random 
walk among traps whose energy depth E > is exponen- 
tially distributed f{E) = e^p{~E/Tg)/Tg where Tg is a 
measure of the disorder. It is easy to show that a — T/Tg 
and A = \T{—a)\a where T is the thermal temperature. 
The goal of this paper is to find the long time behavior 
of {P{x, t)) the probability of finding the particle on x at 
time t averaged over the disorder. For a comprehensive 
mathematical review of the QTM see [l^l ■ 



Q 0.2 0.4 0.6 0.8 1 
{z^} 0.5 0.67 0.81 0.91 0.96 1 

TABLE I: Simple Brownian simulations on a lattice give (z^) 
which according to Eq. fTO]) yield (x^) for the QTM. 

Time in the quenched trap model is t = X]^-oo "i^xTx 
where Ux is the number of visits to lattice point x. Define 



the random variable 77 = t/{SaY^'^ where 

00 

a:— — 00 

When a = 1, 5q is the total number of jumps made 
X]^-oo ~ ^- the opposite limit a — > 0, 5o is the 
distinct number of sites visited by the random walker 
which is called the span of the random walk. We now 
show that in the scaling limit 



the PDF of ?7 is: la,A,i {rj) , 
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where la,A,i{il) is the one sided Levy PDF whose Laplace 
•q ^ u pair is exp(— Au"). Namely the heavy tailed dis- 
tribution of the waiting times Tx determines the statistics 
of 7} through the characteristic exponent a, while the vis- 
itation numbers {ux} provide the scaling through Sa- By 
definition the Laplace rj ^ u transform of the PDF of 77 
is 
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We average with respect to the disorder, namely 
with respect to the independent and identically dis- 
tributed random waiting times Tx, and obtain (e^'"') = 



n; 



of the PDF of waiting times iI}{tx) 
exp(-Aw") - 1 - Au" -I 



where il^iu) is the Laplace transform 
We use i!{u) ~ 



nr=-oo exp 
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Hence the PDF of 77 is a one sided Levy law Eq. ([3]). 
In a longer publication we will complete the proof and 
consider the case where ipiTx) belongs to the domain of 
attraction Levy PDFs (i.e. families of PDFs satisfying 
ip{u) ^ 1 — j4ii" + •••). We now invert the process fixing 
time t to find the PDF of Sa 



nt (Sa) = -{Sa)-^l''-Ha,A,l 
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In the following we explain how to use the operational 
time Sa to obtain the desired diffusion front of the QTM, 
in other words we explain how we get rid of the disorder 
and focus only on Brownian motion. 

Weak subordination breaking. To find {P(x,t)) we fol- 
low six steps: 1. Choose the laboratory time t which 
is a fixed parameter. 2. Use a random number genera- 
tor and draw the stable random variable 77 from the one 
sided Levy PDF la.A,i{il)- With 77 and t determine 
the operational time Sa — {t/v)"- 4. Generate a sim- 
ple symmetric random walk on a lattice (probability 1/2 
for jumping left and right). Stop the process once its Sa 
reaches the operational time set in step 3. 5. Record 




FIG. 1: The PDF Ba{z) exhibits a transition between a 
Gaussian shape when a — >■ 1 to a 1/ shape when a — > 0. 
Simulations of Brownian motion on a lattice yield excellent 
agreement with theoretical predictions Eqs. (|12l llSp without 
fitting. 



the position x of the particle at the end of the previous 
step. 6. Go to step 2. After this loop is repeated many 
times, we generate a histogram of x. The histogram so 
created is identical to {P{x,t)) when t is large. On a 
computer the second step is implemented with a simple 
algorithm provided by Chambers et al 17[. Notice that 
with this exact scheme we have mapped the random walk 
in a random environment to a Brownian motion problem. 
We sec that for quenched disorder the operational time 
is Sa and in this sense subordination is weakly broken: 
the Levy transformation [l^ is still maintained. 

The diffusion front of the QTM. Let Ps^{x) be the 
PDF of X for the simple random walk on a lattice (Brow- 
nian motion) stopped at the operational time Sa- Since 
the QTM dynamics can be separated into two distinct 
processes: Brownian motion with operational time Sa 
(step 4.) and the Levy time transformation (steps 2. 
and 3.) we find 
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where nt {Sa ) is given in Eq. ^ . For the mean field ver- 
sion of the model (i.e. CTRW) replace Sa with the num- 
ber of steps s of the Brownian motion, and then Pg^ (x) 
is Gaussian as well known [l5| . From normal Brownian 
motion we have the scaling behavior x oc (iSq)^/^^^"-*. 
To see this we use: (i) usual Brownian scaling x oc s^/^ 
(ii) 72x within a region |a;| < s^^^ is roughly the num- 
ber of jumps made s divided by the number of sites 
in the explored region Ux oc s/s^/'^ — s^/^. Hence 
Sa oc y/s{nx)°' OC s^^+^^Z^ which gives x oc (5a)^/(^+"). 
This scaling implies 
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FIG. 2: The mean square displacement of the QTM versus 
time. Numerical data matches perfectly the theory based on 
weak subordination breaking (the lines plotted with {z^) in 
Table |T]| and analytical formulas Eq. (|13p for ct = 0.2 and Eq. 
(Unj for a = 0.4,0.6,0.8. 

with Ba{z) a normalized non negative function. Define 
the scaling variable ^ = x/ (t/A-^^") and {P{x,t)) ~ 
9a (C) / {t/A^/°')^ which according to Eq. © is 

poo 

5a (0= / d2/2/TfeS,(e2/Tfr)/„,i,i(y). (9) 

A general relation is found between the moments (|a:;|'') = 
ilZo \x\''P{x,t)dx) of the original QTM and the mo- 
mei^s =/:°^|z|«B„(z)dz 

r ( —2—^ / \ 1^ 

The new content of Eqs. (|9l [TOl) is that once we obtain 
Ba{z) either from theory or simulations of Brownian tra- 
jectories, we have a useful method to obtain exact statis- 
tical properties of the diffusion front. 

Generating Brownian trajectories on a lattice we found 
Ba{z) in Fig. [U which shows an interesting transition 
from a V shape when a — > to a Gaussian shape, which 
we soon analyze analytically. With (z^) given in Table 
H] and Eq. (|10p we get the mean square displacement of 
the QTM (x^). We then favorably compare the predic- 
tions of our theory with simulations of the QTM in Fig. 
[2] (and analytical formulas soon developed). In Fig. [3] 
we show 5a (C) ^-i^d present excellent agreement between 
weak subordination breaking and direct simulation of the 



FIG. 3: The diffusion front of QTM (squares) perfectly 
matches theory based on weak sub-ordination breaking (cir- 
cles) and analytical predictions (lines) Eqs. (|9ll4ll5p . 



QTM. One advantage of our approach is that it is capa- 
ble of dealing with the critical slowing down pointed out 
by Bertin and Bouchaud 16]. Briefly, QTM simulations 
do not converge on reasonable computer time scales for 
say a > 0.8. In contrast weak subordination breaking 
scheme quickly converges since it is based on Brownian 
motion and there is no need to generate disordered sys- 
tems. More importantly we now analyze Brownian mo- 
tion analytically, obtain Ba{z) in two important limits 
and then with Eqs. (HllTni) provide solutions to the QTM. 

The limit a — corresponds to strong disorder. To 
find Bq(z) we consider Brownian motion stopped at 
"time" So where as mentioned Sq is the span of the ran- 
dom walk. Consider Psg {Sq — n) where x = Sq — n > 
and for simplicity we start with n = 1. The Brownian 
particle after the first step can be either on a; = 1 or 
X — —1. If it is on a; = —1 it must travel a distance 
Sq to reach its destination x = Sq ~ 1 and the span is 
S{). On the other hand if it jumps to a; = 1 the distance 
the particle must travel is 5o — 2 and the span must still 
be Sq. Hence Ps„ {Sq - 1) = [Ps, {Sq) + Ps„ {Sq - 2)]/2. 
More generally 

Ps„ {Sq 'n) = ^ [Ps„ (5o - n - 1) + Ps„ (5o - n -f 1)] , 

(11) 

and for the boundary term Ps^-,{Sq) = 
[Pso{Sq-1) + Pso-i{Sq-1)]/2. Eq. ^ is eas- 
ily solved Pso {x) = for —Sq < x < Sq and 
X G Z. In the limit 5o ^ 1 we have for the scaled 
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variable z = xjSf^ the V shape PDF (see Fig. [T]) 
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This V shape reflects the tendency of a Brownian particle 
to reach a large span Sq when it is far from the origin. 

According to Eq. pO)) the even moments (x^') for 
the random walk in the QTM are given once we obtain 
{z^i). In the limit a we find using Eq. \Y1\ {z^'') — 
2 z^izdz 
small a 



(1 + g) ^ hence with Eq. ([T0| we have for 
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which is tested in Fig. [21 Inserting (z^"^) = (1 + g)~^ 
in Eq. dTUl) we obtain the moments {x^'^) of the QTM. 
Straight forward analysis then gives 
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where Ei[^) = J^{e~^/t)dt is the tabulated exponential 
integral. This scaling function was obtained by C. Mon- 
thus 14] using an RG method which is exact in the limit 

Approaching the weak disorder limit a — ^ 1. For 
a = 1 we have Si = J2T=-oo ~ namely Si 
is non random since it is equal to the number of 
steps made. Therefore when a is close enough to 1 
we may neglect fluctuations and Sa = {Sa}- In a 
longer publication we will show that (Sa) = CaS^~ . 
and Cc. = 2("+3)/2p(j_^Q,/2)/[0f(i + Q,)] [hints: find 
{{nx)") with the first passage time PDF for Brown- 
ian motion and renewal theory, for the limit a — >■ 
consult fisll]. As well known the PDF of finding the 
Brownian particle on x after s jumps is the Gaussian 
Ps{x) — exp (— a;^/2s) /\/2tts hence change of variables 
to Sa gives 



exp 



Baiz) 



2 



27r/(C„) 



1/2 



(15) 



It follows that (z^) ^ (Ca) hence for the QTM 
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In Fig. [T]_Bq,(z) obtained from Brownian simulations is 
favorably compared with Eq. (1151) for a = 0.9. Surpris- 
ingly, as we show in Fig. [2] Eq. (|16p works very well 



even for a = 0.4. With Eqs. (jH [TS]) and steepest descent 
method we find for ^ >> 1 
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which approach the expected normal Gaussian limit 
when a — >■ 1. In the opposite limit ^ << 1 (?a (0 ^ 
l/V2^-2("-i)/2 [(1 -I- a)/a] {Ca/T [(1 - «)/2]}r + ' ' ' ■ 

The method presented here is not limited to one dimen- 
sion neither to unbiased motion. For example consider 
the QTM on a regular lattice in three dimensions. From 
ordinary Brownian motion we expect Sa ^ CaS, where 
now Sa is non random and the (non trivial) parameter 
Ca will depend on the lattice structure. In principle once 
Ca is determined one can map the QTM with an equation 
like Eq. ([7]) to an ordinary Brownian motion. This will 
yield mean field CTRW dynamics which is still non triv- 
ial since the transformation depends on the parameter 
Ca- Note that the QTM dynamics describes also certain 
models of random walks on random geometries, e.g. frac- 
tal comb structures and naturally our approach captures 
also these cases. 

Finally, weak subordination breaking scheme is not 
only a new approach which deals with anomalous diffu- 
sion in systems with quenched disorder. Our method can 
be used to solve other aspects of dynamics in disordered 
systems like aging and weak ergodicity breaking. 
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